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The anomaly in the low energy distribution of quasi-elastic neutrino events reported by the Mini- 
BooNE collaboration is discussed. We show that the observed excess of electron-like events could 
originate from the production and decay of a heavy neutrino (uh) in the MiniBooNE detector. The 
Ph is created by mixing in v^, neutral-current interactions and decays radiatively into vy due to a 
transition magnetic moment between the Uh and a light neutrino u. The energy measured in the 
detector arises from the subsequent conversion of the decay photon into a e'^e~ pair within the 
detector volume. The analysis of the energy and angular distributions of the excess events suggests 
that the Vh has a mass around 500 MeV and the lifetime r^^ < 10"^ s. Existing experimental data 
are found to be consistent with a mixing strength between the Vh and the Vp, of |?7,jh P — (1 — 4) x 10~^ 
and a Vh transition magnetic moment of ^tr — (1 — 6) x 10~^/iB. Finally, we discuss the reason why 
no significant excess of low energy events has been observed in the recent antineutrino data. 

PACS numbers: 14.80.-j, 12.60.-i, 13.20.Cz, 13.35.Hb 



The MiniBooNE collaboration, which was searching 
for the LSND signal from neutrino oscillations at FNAL, 
has observed an excess of low energy electron-like events 
in charge-current quasi-elastic {C'CQE) neutrino events 
over the expected standard neutrino interactions [^. 
This anomaly has been recently confirmed by the finding 
of more excess events While the collalDoration has 
not yet clarified the origin of the excess, several mod- 
els involving new physics were considered to explain the 
discrepancy 

In this work we show that the excess could be explained 
by the production and decay of a heavy neutrino (vh)- 
Such type of neutrinos are present in many interesting 
extensions of the Standard Model, such as GUT, Super- 
string inspired models, Left-Right Symmetric models and 
others. The massive neutrino decays were also considered 
to explain the LSND signal [3]. 

The neutrino weak flavor eigenstates (I'e, v^, i^t,---) 
can be different from the mass eigenstates 
(j/i, 1^2, 1^3, J^4...), but they are related to them, 
in general, through a unitary transformation. A 
generalized mixing: 

vi = '^UiiUi; Z = e,^,T, i = 1,2,3,4, ... (1) 

i 

results in neutrino oscillations when the mass differences 
are small, and in neutrino decays when the mass differ- 
ences are large. If the Vh exists, it could be a component 
of f^, and, as follows from Eq.([T]), it would be produced 
by any source of according to the proper mixing | U^h\'^ 
and kinematic constraints. The Vh can be Dirac or Ma- 
jorana type 0] and can decay radiatively into v^, if there 
is a non-zero transition magnetic moment between the Vh 
and a light neutrino Vi [6]. 

The MiniBooNE detector is described in details in Ref. 
. It is using an almost pure beam originated from 
the tt"*" decays in flight, which are generated by 8 GeV 
protons from the FNAL booster. The detector consists 
of a target, which is a 12.2 m diameter sphere filled with 



800 t of mineral oil, surrounded by an outer veto region. 
The Cherenkov light rings generated by muon, electron 
and converted photon tracks are used for the reconstruc- 
tion of the events. The resolutions reached on the vertex 
position, the outgoing particle direction and the visible 
energy are 20 cm, 4 degrees and 12%, respectively for 
CCQE electrons [sj]. The v^^ beam is peaked around 
~ 600 MeV, has a mean energyof ~ 800 MeV and a 
high energy tail up to 3 GeV @- 

An excess of AiV =128.8±20.4 ± 38.3 electron-like 
events has been observed in the data accumulated with 
6.64 X 10^° protons on target. For the following discus- 
sion several distinctive features of the excess events are 
of importance 2]: a) the excess is observed for single 
track events, originating either from an electron, or from 
a photon converted into a e^e~ pair with a typical open- 
ing angle ~ me/E^+^- < 1 degree (for E^+^- > 200 
MeV), which is too small to be resolved into two separate 
Cherenkov rings (here, TOe,£^e+e- ^'I'e the electron mass 
and the e"'"e~ pair energy); b) the reconstructed neutrino 
energy is in the range 200 < E^^ < 475 MeV, while 
there is no significant excess for the region E^^ > 475 
MeV. The variable E^^ is calculated under the assump- 
tion that the observed electron track originates from a 
i/g QE interaction; c) the visible energy E^is is in the 
narrow region 200 < Eyis < 400 MeV for events with 
E^^ > 200 MeV; d) the angular distribution of the ex- 
cess events with respect to the incident neutrino direc- 
tion is wide and consistent with the shape expected from 
Vf,CC interactions. To satisfy the criteria a)-d), we pro- 
pose that the excess events are originated from the de- 
cay of a heavy neutrino Vh- The v^s are produced by 
mixing in t/^ neutral-current (NC) QE interactions and 
deposit their energy via the visible decay mode i^J, 
as shown in FiglU with the subsequent conversion of the 
decay photon into e+e^ pair in the MiniBooNE target. 
To make a quantitative estimate, we performed simpli- 
fied simulations of the production and decay processes 
shown in Fig[TJ In these simulations we used a energy 
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FIG. 1: Schematic illustration of the NCQE production and 
decay of heavy neutrino. 



spectrum parametrized from the reconstructed v^CCQE 
events @. Since in the MiniBooNE experiment VhS de- 
cay over an average distance of < 5 m from the produc- 
tion vertex, the sensitivity is restricted to the mass range 
100 - 600 MeV and to vh hfetimcs t^^ < 10"'^ s. 

In Fig[21ll]the distributions of the kinematic variables 
E^^,Eyis and cos{Qj^) for the i^h — > 71^ events are 
shown for m^^ = 400 and 600 MeV and r^,^ = 3 x 10"® 
and lO^^^s. These distributions were obtained assuming 
that the e'^e~ pair from the converted photon is mis- 
reconstructed as a single track from the v^QE reaction. 

Simulations are in reasonable agreement with the ex- 
perimental distributions. For instance, for the distribu- 
tion shown in FigEl the comparison with MiniBooNE 
data yields a of 10.2 (17.2) for 8 DF corresponding 
to 27% (~ 5%) CL. for m^,^ = 400 (600) MeV and 
Tuf^ = 3 X 10^^ s. The simulated excess events, shown 
in Fig (HI are mainly distributed in the narrow region 
200 < Ey,s < 400 MeV. The fraction of events in the 
region 200 < E^is < 400 MeV is ~ 70%. The remaining 
events are distributed over the region 400 < Eyis < 1200 
MeV, where they can be hidden by the low statistics. 
The simulations showed that the shape of the E^^ and 
Eyis distributions is sensitive to the choice of the mass 
and lifetime: the shorter the Vh lifetime the broader the 
visible energy spectrum. The best fit results suggest that 
the i/fi mass is in the region 200 < m^^ < 600 MeV and 
the lifetime is r^^^ < 10~^ s. The estimate of the mixing 
parameter IC/^^P was performed by using the following 
relations. For a given flux of heavy neutrinos, <I>(j/^), the 
expected number of the decays in the MiniBooNE detec- 
tor is given by AN = / ^{iyh)PdecPconvedE^^dV , where 
Pdec and Pconv are the probabilities of the Vh decay and 
the photon conversion in the detector, e is the overall 
detection efficiency, and the integral is taken over the 
detector fiducial volume. 

The flux ^{vh) was estimated from the expected num- 
ber of the Vfj^NC events times the mixing iJ/^^/iP, taking 
into account the threshold effect due to the heavy neu- 
trino mass. The total number of reconstructed v^j^CC 
events in the detector [9;] was used for normalization. 
The probability of the heavy neutrino to decay radia- 
tively in the fiducial volume at a distance r from the pri- 
mary vertex is given by Pdec = [l-exp( )] ^^"r'^^''^ 
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FIG. 2: Distributions of the excess events from the Vh 'yv 
decay reconstructed as UeCC events as a function of E^^ for 
\U^,hf = 1.5 X 10"^ and a) = 400 and r^^ = 3 x 10"* s 
{fitr = 2 X 10~"Vs) (solid); b) = 400 and r^^ = 10~^° s 
(/itr = 3xl0~^/is) (dashed); c) m^^ = 600 and Ti,^ = 3xl0~* 
s (dotted). The dots are experimental points for the excess 
events in the MiniBooNE detector. Error bars include both 
statistical and systematic errors 0]. The comparison of the 
distributions with the experimental data yields a of 10.2, 
11.2, and 17.2 for 8 DF corresponding to 27%, 24%, and ~ 
5% C.L. for a),b) and c), respectively. 



where the last term is the branching fraction Br{yh — > 
^v) ~ 1 (see below). Taking into account the ratio 
v^,NCQE/v^,CCQE ~ 0.43, the number of v^^CCQE 
events observed P, |3| and assuming that almost all 
Vh — > IV decays occur inside the fiducial volume of the 
detector, we estimate the iJ/^^/iP to be in the range 

|C/^^P ~ (1-4) X 10-3. (2) 

This result is mainly defined by the uncertainty on the 
number of excess events. Eq.Q is valid for the mass 
region 400 < to^^ < 600 MeV. The lower limit is set to 
400 MeV to avoid stringent constraints on |?7^/ip for the 
mass region ^ 400 MeV from experiments searching 
for a peak from n,K — > ^ + decays [lo| . The Vh 
lifetime due to a transition moment fitr is given by 

and for t-^i, < 10"'' s results in fitr > 10~^°iib- The to- 
tal Vfi decay width is Ttot = ^{i^h i^"f) + Sri, where 
r(i^?i ^ i^"f) is the Vh 7*^ decay rate, and EF^ is the 
sum over decay modes whose decay rate is proportional 
to the square of the mixing IC/^^P- The dominant con- 
tribution to SFi comes from Vh —>■ v^ee,v^T:^ ^jltt 
decay;s, for which the rate calculations can be found, e.g. 
in 14 For m^^ ~ 500 MeV and /it,. > 10~i°/xb, we 
found that the radiative decay is dominant, Br(vh — > 
71^) > 0.5. For example, for fitr = 10~^/i_B and 
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FIG. 3: Distributions of the excess events from the Vh — > 'yv 
decay reconstructed as Ue.CC events as a function of Evia for 
B?^ > 200 MeV. The comparison of the distributions with 
the experimental data yield a of 9.7, 10.3, and 16.8 for 8 
DF corresponding to 28%, 27%, and ~ 5% C.L. for a), b) and 
c), respectively. The legend is the same as in Fig|2] 



FIG. 4: Distribution of the excess events from the Vh decay 
reconstructed as Ve.CC events as a function of cos 671^ for 
300 < < 400 MeV. The comparison of the distributions 

with the experimental data yields a of 11.6, 11.1, and 15.6 
for 8 DF corresponding to 23%, 24%, and ~5% C.L. for a), b) 
and c), respectively. The notations are the same as in Fig[2] 



W^h]"^ = 1.5 X 10 the expected ratio of decay rates 
for -fv : ^JLTT : e^v : fifiiy is 0.984 : 0.011 : 0.0016 : 
0.00067. 

One may wonder if the mixing strength of Eq.Q is 
consistent with the results of previous searches for 
decays. The i/h mass region around 500 MeV was cov- 
ered by many experiments p^ . [ill , [l^ . However, none 
of these experiments has reported a bound on the mix- 
ing strength or on the combination \Uf_ih\'^ f^tr, for 
the radiative Vh neutrino decay. The best Umit 
\Uiih\'^ ^ 10^^ for the mass region m^f^ ~ 500 MeV 
was derived from a search for Uh ^tt, ujxv, nev de- 
cays in the NuTeV beam dump experiment [l3| (see also 
[3 d, It was assumed that these decay modes 
are dominant and that the Vh is a relatively long lived 
particle, i.e ^^'"'^ <^ 1, where L ~ 1.4 x 10^ m is the 

P''h '''"h 

distance between the target and the detector. Consider 
now our case with ji/^hp = (1 — 4) x 10~^, m^^ = 500 
MeV and Utr — 10~^/iB. This gives the Vh lifetime 
T^^ — (1.5 — 1.4) X 10^^ s. Due to the larger mix- 
ing the i^h flux at the target would increase by a factor 
~ (1 — 4) X 10'^. However, taking into account the attenu- 
ation of the flux due to the rapid decay of v^s, the total 
number of signal events in NuTeV would decrease by a 
factor (10—3) compared to the number of events expected 
for a long lived Vh's produced and decaying through the 
mixing |t/;i/ip = 10~^. In this estimate the average Vh 
momentum is < p^^, >~ 100 GeV and the decay region 
length is ^ = 34 m 13]. Finally, we find that for 

Mtr > IO^Vb (4) 

the NuTeV limit is not constraining mixing of Eq.(l2]). 
Note that a short Vh lifetime is also necessary to avoid the 



constraints coming from cosmological and astrophysical 
considerations p^ . 

The best limit for the large mixing (short lifetimes), 
C/p/iP < 10"'^ for the masses around 500 MeV, was de- 
rived by CHARM-II from a search for the Vh production 
and — > iifiv decays within their detector [l5l |. Tak- 
ing into account the Vh — * 7^^ decay, the < >~ 24 
GeV and Z = 35 m, we find that the number of the 
expected fifxiy signal events in CHARM-II is iV^^^ ~ 

0.46 X 10"^^ (J^^^T^. For /ij^ > 3 x lO'Ve and for mix- 
ing of Eq.jSI) this results in N^^,^ < 0.05 - 0.8 events. To 
evade the CHARM-II limit for the region 10~^ < Utr < 
3 X 10^^/is the mixing is required to be in a slightly more 
restricted range 10^"^ < \Uf^h\'^ < 1.5 x 10^^ compare 
to that of Eq.Q. For example, for utr = 2 x lO^^/is 
the allowed range is 10~^ < \Uf,h\'^ < 3 x 10"^. Thus, 
we see that most of the allowed {fitr', \ U^h\'^) parameter 
space corresponding to Eas. (l2l4p is not constrained by 
the CHARM-II limit. 

Consider now bounds from LEP experiments . For 
the mass region around 500 MeV, the model independent 
limit from the searches for the Z vvh decay is |C/phP ^ 
10^^, (see e.g. [11]) which is compatible with Eq.®. 
Direct searches for radiative decays of an excited neutrino 
V* — !■ 71/ produced in Z ^ v*u decays have been also 
performed |l3]. The best Hmit from ALEPH is Br{Z 
vv*)Br{v* ^v) < 2.7x 10~^ [11]. As the experimental 
signatures for the v* — > and Vh — >■ 7^^ decays are the 
same, we will use this bound for comparison. The number 
of expected Vh — > 7^^ events in ALEPH is proportional to 
Br{Z ^ vvh)Br{vh ^ 7'^)[1 -exp(--5^)], with I ~ 1 
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m and p,^,^ ~ 45 GeV. Taking into account Br(z^vv) — 
\Uf_,h\^ and using Eq.dS]), we find 

IC^mI'x (—)'< 3.5x10-20. (5) 

Using Eq.([21) results in /ij,. ^ (6 — 3) x IQ-^/is, which is 
consistent with Eq.Q. 

The hmit on the ^tr between the and the v^^ has 
been obtained in Ref . [20| . based on the idea of the Pri- 
makoff conversion v^Z ^ UhZ oi the muon neutrino into 
a heavy neutrino in the external Coulomb field of a nu- 
cleus Z, with the subsequent Vh — > 71/ decay. By us- 
ing the results from the NOMAD experiment [2ll . [23 |. 
a model-independent bound ^^l' < 10~®^b was set for 
the Vh masses around 500 MeV (see Table 1 and Fig. 2 in 
Ref.[20j). which is also consistent with Eq.([4]). 

The low statistics anti-neutrino {y^) data collecte d by 
the MiniBooNe seem to show no low-energy excess |23| . 
An analysis of these data within the framework discussed 
above suggests that the excess is not seen due to the lower 
77^ energy. Indeed, the 77^ flux peaks at ^ 400 MeV and 
has a mean energy of ^ 600 MeV [§]■ If the Vh mass is 
around 500 MeV, the Vh production is kinematically sup- 
pressed for energies below the mean energy. Instead of 
the expected excess of ~ 40 events [23:] , a smaller excess 
of ~ 23 events is expected in the antineutrino data. 

In summary, we see that the interpretation of the Mini- 



BooNe anomaly based on the production and visible de- 
cay of a heavy neutrino is compatible with all the four 
constraints a)-d). The shape of the excess events in sev- 
eral kinematic variables is found to be consistent with the 
distributions obtained within this interpretation. The 
reason why the excess is not observed in the recent an- 
tineutrino data is clarified. A definite conclusion on 
the presence of Vh — > 17^7 events can be drawn when the 
Vfj^ statistics is substantially increased. Our results for 
the mixing strength |C/p/iP ~ (1 — 4) x lO^^^ and for the 
magnetic moment ^tr — (1 ~ 6) x lO^^/is are compat- 
ible with the results from previous experiments. Values 
of y,tr larger than 10~^°/^b could be obtained e.g. in the 
framework of the Zee model . Our analysis gives a cor- 
rect order of magnitude for the parameters P and ^tr 
and may be improved by more accurate and detailed sim- 
ulations of the MiniBooNE detector, which are beyond 
the scope of this work. We note that an analysis of the 
excess of events due to the Vh — > 'yv decay may also be 
possible with existing neutrino data. New results could 
be obtained with NOMAD [llLSciBooNE [H and K2K 
near detectors [1^, see also[2^|. The author thanks S. 
Brice, A.D. Dolgov, S.H. Hansen, D.S. Gorbunov, N.V. 
Krasnikov, V.A. Matveev, V.A. Rubakov, M.E. Shaposh- 
nikov, and R. Van de Water for useful discussions and/or 
comments, and M. Kirsanov, R. Petti and D. Sillou for 
help. This work was supported by Grant RFBR 08-02- 
91007-CERN. 
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